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I . SUMMARY 

During  the  past  several  years.  Systems,  Science  and 
Software  (S3)  personnel  have  been  actively  engaged  in  a com- 
prehensive program  involving  computer  modeling  of  the  non- 
lihear  processes  that  characterize  underground  nuclear  ex- 
plosions, propagation  of  the  resultant  3tress  waves  through 
realistic  earth  structures  and  prediction  of  the  ground 
motion  recorded  at  teleseismic  distances  from  an  explosion. 
The  objectives  of  the  subject  project  are  to  employ  these 
modeling  and  predictive  capabilities  in  a systematic  examina- 
tion of  the  effects  of  variations  in  source  and  emplacement 
parameters  on  seismic  signals  from  underground  explosions, 
and  to  investigate  methods  for  utilizing  the  general  charac- 
teristics of  seismic  waveforms  to  obtain  reliable  yield  esti- 
mates for  explosions. 

The  technical  phases  necessary  to  accomplish  the  ob- 
jectives of  this  project  are  as  follows: 

1.  Conduct  a systematic  theoretical  examination 
of  material,  source  and  emplacement  parameters 
which  affect  yield-magnitude  relationships 
and  compare  the  theoretical  predictions  to 
actual  ooservations. 

2.  Determine  and  express  uncertainties  cZ  yield 
estimates  in  terms  of  uncertainties  in  gross 
earth  structure,  near  source  material  proper- 
ties, and  local  source  and  receiver  structure. 

Major  accomplishments  during  the  first  three  months  of 
this  project  were  realized  in  several  different  areas  of  re- 
search. In  particular,  a very  comprehensive  investigation 
of  the  sensitivity  of  seismic  magnitudes,  which  are  used  to 
estimate  explosion  yields,  to  variations  in  near  source 
material  properties  was  carried  out.  By  means  of  this  study 
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several  parameters  were  isolated  and  systematically  varied 
to  determine  their  relative  influence  on  seismic  magnitudes. 

A discussion  and  summary  of  the  results  of  this  study  form  a 
subsequent  part  of  this  report.  A technical  report  describ- 
ing the  results  of  this  parameter  study  in  greater  detail  is 
in  preparation  (Wray,  et  al . ^ ) . 

Another  area  of  particular  importance  is  the  exercise 

of  our  ground  motion  predictive  capabilities  for  specific 

NTS  underground  explosions.  The  near  field  coupling  and  the 

[21 

equivalent  elastic  source  (Cherry,  et  al . J)  were  computed 
for  the  recent  explosion,  Mast,  detonated  at  Pahute  Mesa. 
Information  on  certain  material  properties  appropriate  for 
the  shot  medium  were  provided  in  a CEP  document.  Synthetic 
seismograms  will  be  computed  at  receiver  sites  specified  by 
the  Project  Officer.  There  are  several  important  questions 
that  we  will  try  to  answer  as  a result  of  this  experiment. 

For  instance,  what  pre-shot  measurements  of  the  near  source 
medium  are  required  in  order  to  be  able  to  predict  the  seismi- 
cally  determined  yield  to  within  some  specified  range?  The 
results  of  this  experiment  have  important  implications  for  a 
Limited  Threshold  Test  Ban  Treaty. 

Considerable  effort  was  spent  on  code  development 
during  the  past  three  months.  A new  tension  fracture  model 
was  developed  in  the  shock  code  used  for  modeling  the  non- 
linear regime  of  a nuclear  explosion.  Additional  calcula- 
tional  tools,  one  that  will  facilitate  modeling  of  two- 
dimensional  phenomena  such  as  cratering,  and  another  that 
will  provide  a more  accurate  generalized  ray  theory  approach 
to  computing  amplitudes  of  seismic  waves  traveling  through 
the  earth's  upper  mantle,  were  added  to  the  S3  library  of 
computer  codes. 
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II.  INTRODUCTION 


As  stated  in  the  previous  section  the  objectives  of 
the  subject  research  project  are  to  utilize  existing  computa- 
tional capabilities  to  examine  the  effects  of  various  source 
and  emplacement  parameters  on  seismic  signals  from  under- 
ground explosions,  and  to  devise  and  evaluate  methods  for 
utilizing  the  general  characteristics  of  seismic  waveforms  to 
obtain  reliable  estimates  of  explosion  yields.  In  order  to 
realize  these  objectives,  activity  on  this  program  during  the 
first  three  months  of  the  contract  period  was  concentrated  in 
the  following  areas: 

1.  Study  of  the  dependence  of  teleseismic  magnitude 
on  near  source  material  properties. 

2.  Prediction  of  teleseismic  ground  motion  for 
specific  NTS  explosions. 

3.  Improvements  in  computational  capabilities. 

The  plan  of  the  remainder  of  this  report  is  a technical  discus 
sion  of  each  of  these  three  research  areas,  followed  by  a 
section  summarizing  the  most  important  results  obtained  to 
date . 
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III.  TECHNICAL  DISCUSSION 


3.1  INTRODUCTION 

A key  requirement  for  predicting  both  body  and  surface 
wave  teleseismic  ground  motion  from  a nuclear  explosion  is 
the  specification  of  the  equivalent  elastic  source  of  the  ex- 
plosion. Computational  techniques  and  material  models  have 
been  developed  at  S3  which  permit  a prediction  of  the  equiva- 
lent elastic  source  given  the  yield  of  the  explosion  and  the 
material  properties  of  the  near  source  geologic  environment. 
In  the  case  of  a spherically  symmetric  explosion  source  a 
convenient  representation  of  the  equivalent  elastic  source 
is  given  in  terms  of  the  steady  state  value  of  the ^reduced 
displacement  potential  (RDP)  , ^ (Cherry,  et  al.  ). 

A particularly  important  relationship  between  the  tele 
seismically  determined  bodywave  magnitude,  mb,  and  ^ is 


m,  ~ log 

D 


(°°) 

c 


(3.1) 


[2] 


where  C is  the  near  source  sound  speed  (Cherry,  et  al. 
This  relationship  will  be  used  often  in  the  following  sub- 
sections  of  this  report. 


3.2  dependence  OF  teleseismic  magnitude  on  near  source 

MATERIAL  PROPERTIES 

The  basic  objective  of  the  study  presented  in  this 
section  is  to  determine  the  sensitivity  of  seismic  yield- 
magnitude  relationships  to  parameters  related  to  the  explosive 
source  and  the  near  source  material  properties.  Two  distinct 
source  descriptions  have  been  employed  and  the  differences  are 
discussed  in  subsection  3.2.2.  The  remaining  parameters  are 
numerous  and  are  all  related  to  material  description.  Since 
our  constitutive  models  are  fairly  complex,  great  care  has 
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been  exercised  to  conduct  the  sensitivity  study  in  a syste- 
■i  matic  and  meaningful  manner.  A general  discussion  of  the 

organization  of  the  sensitivity  study  and  the  meaning  of  the 
various  independent  material  parameters  is  presented  in  the 
following  paragraphs  of  this  subsection. 

Figure  3.1  depicts,  in  flow  chart  form,  the  computa- 
tional organization  followed  in  the  one-dimensional  para- 
meter study.  The  study  begins  with  a material  composed  of  a 
mixture  of  tuff  and  water  which  may  contain  air-filled  voids. 

The  properties  of  the  water  and  compacted  tuff  constituents 
are  held  constant  for  all  calculations.  The  water  is  described 
by  an  equation  of  state  formulation  developed  by  Bjork^and 
Gittings[3]  and  the  formulation  given  by  Riney,  et  al^  is 
used  for  tuff  with  a grain  density  of  2.4  gm/cm.  The  mass 
fraction  of  water  in  the  mixture,  fw»  is  a natural  parameter 
to  select  for  initial  variation  and  this  choice  is  indicated 
in  Fig.  3.1  which  contains  a "flow  chart"  of  the  sensitivity 
study.  The  water  content  study  was  conducted  by  varying  fw 
in  a material  with  no  air-filled  voids;  the  procedure  and  re- 
sults are  discussed  in  subsection  3.2.1. 

The  remaining  material  parameters  are  varied  one  at  a 

time  while  holding  the  water  mass  fraction  constant  at  0.17. 

Four  of  these  material  parameters  are  directly  related  to  the 

[ 5 I 

porosity  model  which  is  discussed  by  Cherry,  et  al.: 

(1)  the  air-filled  volume  fraction,  4>  o ; (2)  the  elastic  limit 

of  the  pressure,  Pe;  (3)  the  crush  or  compaction  pressure,  Pc? 
(4)  the  bulk  modulus,  k.  It  should  be  recognized  that  k may 
not  be  varied  with  complete  freedom,  after  the  constituent 
properties  and  the  mass  fraction  of  water  are  fixed,  without 
violating  consistency  requirements.  For  the  second  indepen- 
dent elastic  property  of  the  porous  mixture  we  select  the 
shear  or  rigidity  modulus,  y.  The  compressional  wave  velocity, 

c,  equals  /k  + 4/3  I[ 

P 
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Three  additional  material  parameters  are  contained  in 
the  pressure  and  energy  dependent  model  for  the  yield  strength. 
The  yield  strength,  Y,  is  given  by 


Y = 


P / P 

Vo  + Y*  5-  l2  ■ E 


m 


)|  i1  • y 


, p < p 


■ (Y. + Y*)  i1  - k) 


, , p > p 
e / - m 

m 


m 


e < c 


m 


(3.2) 


= 0,  e > e 
' — m 


where  Y is  the  zero  pressure  value  of  the  yield  strength, 
o 

YMAX  = + Ym  the  max^mum  yield  strength,  P is  the  pres- 

sure, Pm  is  the  pressure  at  maximum  yield  strength,  e is  the 
energy  density  and  em  is  the  energy  density  at  incipient 
melt.  The  melt  energy  is  associated  with  the  tuff  and  is 

therefore  not  varied  leaving  the  independent  parameters  Y , 

0 


ymax'  and  P 


m 

One  last  parameter,  the  overburden  pressure,  P , is 

o 

included  in  the  one-dimensional  sensitivity  study.  Even  though 
the  overburden  pressure  is  a depth  of  burial  parameter  it  is 
included  as  part  of  the  parameter  study. 


3.2.1  Water  Content  Parameter  Study 

The  effects  of  water  content  on  seismic  coupling  were 
studied  using  SKIPPER, ^ a one-dimensional  Lagrangian  non- 
linear wave  propagation  code.  The  reduced  displacement  poten- 
tial (RDP)  was  calculated  for  five  fully  saturated  tuffs 
having  different  mass  fractions  of  water  (f^) . Fully  saturated 
tuffs  were  chosen  for  this  study  in  order  to  eliminate  the  ef- 
fects of  air-filled  voids. 

The  tuffs  were  modeled  using  the  Tabular  Array  of  Mix- 
tures Equation  Of  State  (TAMEOS) . For  a given  mass  fraction 
of  water,  the  TAMEOS  scheme  mixes  dry  rock,  modeled  by  a Mie- 
Gruneisen  equation  of  state,  with  water,  assuming  that  the 
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components  of  the  mixture  are  in  pressure  equilibrium.  Thus, 
the  pressure,  specific  volume,  and  specific  internal  energy 
of  the  mixture  <Pm,  V , and  e^  respectively)  are  given  by  the 
following  equations: 


P = P = P 
m r w 


V = (1  - f )V  + f V 
m w r w w 


e = (1  - f ) e + f e 
m v w r w w 


(3.3) 


where  the  subscripts  r and  w refer  to  the  rock  and  water 
components. 

Five  calculations  were  made  for  tuffs  with  material 

properties  as  given  in  Table  3.1  and  for  a nominal  device  yield 

of  20  tons.  Adding  water  decreases  both  the  density,  p^,  and 

bulk  modulus  Kw  of  the  mixture.  Since  the  shear  modulus  was 
M 

held  constant  at  40  kbars  for  all  calculations,  the  variation 

in  longitudinal  sound  speed  Cfl  reflects  only  the  change  in 

bulk  modulus.  Figure  3.2  represents  graphically  the  results 

shown  in  Table  3.1.  RDP  decreases  as  more  water  is  added  into 

the  mixture.  RDP/C  , a measure  of  the  seismic  amplitude,  is 

o 

also  plotted.  This  also  decreases  as  water  content  is  in- 
creased. 

The  decrease  in  seismic  coupling,  as  the  amount  of 
water  in  the  water-rock  mixture  is  increased,  is  to  be  ex- 
pected. Since  water  is  more  compressible  than  rock,  the 
amount  of  energy  dissipated  in  the  form  of  heat  via  shock 
loading  increases  with  increasing  water  content,  resulting  in 
less  energy  available  for  conversion  into  seismic  radiation. 

The  column  headed  Am  in  Table  3.1  gives  an  indica- 
tion of  the  change  in  teleseismic  magnitude  as  a function  of 

water  content.  The  calculation  for  f = 0.19  was  chosen  as 

w 

the  reference  point.  The  calculation  of  Am  values  is  based 
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on  Eq.  (3.1)  and  is  given  by 


Am  = m - m = log 
1 2 


<P  (00) /C 

1 0 1 

v rt/c 

2 0 2 


where  4>  (°°)  and  ^ (°°)  correspond  to  two  sets  of  maLerial 
1 2 

properties  in  Table  3.1  and  likewise  for  the  sound  speeds  C 

and  C . The  significant  point  about  the  Am  values  in 
2 

Table  3.1  is  that  for  a variation  of  the  mass  fraction  of 
water  of  19  percent  the  change  in  magnitude  is  only  0.09  units. 
As  we  shall  see  in  the  next  subsection  of  this  report,  there 
are  several  other  parameters  that  exert  a stronger  influence 
on  magnitude. 


3.2.2  The  Effect  of  Source  Description,  Unconfined  Yield 
Strength,  Seismic  Velocity  and  Shear  Modulus  on 
Seismic  Magnitudes  ~ 

The  effect  of  source  description  and  material  properties 

on  the  RDP  were  determined  by  performing  calculations  with  the 

[*21 

SKIPPER  code.  A recent  S3  report  J identifies  and  describes 
a number  of  material  parameters  that  strongly  influence  tele- 
seismically  recorded  ground  motion  from  an  underground  explo- 
sion and  presents  the  results  of  a sensitivity  study  performed 
on  those  parameters.  The  results  of  that  study  are  presented 
here  in  Table  3.2  for  convenience.  All  calculations  in  the 
table  were  performed  with  a near  source  sound  speed,  C^,  of 
2.4  km/sec,  a shear  modulus,  y,  of  40  kbar,  an  unconfined  yield 

strength,  Y , of  zero,  a water  mass  fraction,  f , of  0.17,  and 
o w 

an  explosive  yield  of  20  tons. 

, [2] 

The  source  description  employed  in  the  previous  study 
made  use  of  a subroutine  which  computed  the  radial  stress  on 
the  cavity  boundary  via  an  isentropic  relation  for  a gamma  law 
gas.  In  order  to  avoid  excessively  high  pressures  in  the 
earth  material  the  maximum  pressure  obtained  for  the  cavity 
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SUMMARY  OF  PARAMETER  STUDY  USING  ORIGINAL  SOURCE  DESCRIPTION 
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was  not  applied  to  the  boundary  instantaneously.  Instead 
the  cavity  pressure  was  required  to  approach  the  maximum 
value  in  small  increments  taken  over  approximately  50  cycles. 

For  the  new  sensitivity  calculations  reported  herein, 
a more  direct  procedure  is  employed  for  the  source  descrip- 
tion. The  total  initial  energy  of  the  cavity  is  determined 
directly  from  the  device  yield.  The  initial  cavity  radius 
(radius  of  vaporization)  varies  as  the  one-third  power  of  the 
energy  yield  and  is  determined  by  the  methods  discussed  by 
Butkovich. ^ The  initial  energy  density  of  the  cavity  is 
of  course  determined  by  the  total  energy  divided  by  the  total 
mass.  The  cavity  material  is  described  by  a Tillotson  equa- 
tion of  state  which  limits  to  an  ideal  gas  law  at  high  ener- 
gies. A single  spherical  zone  is  used  to  model  the  cavity 
region;  the  cavity  zone  is  fully  coupled  to  the  adjacent  un- 
vaporized earth  material  zones  via  the  conservation  laws  of 
hydrodynamics.  No  ramping  procedure  is  employed.  Several 
calculations  have  shown  that  this  procedure  conserves  energy 
from  the  source  to  within  a few  percent  and  that  more  energy 
is  coupled  into  the  ground  than  with  the  old  source  descrip- 
tion; calculated  seismic  magnitudes  are  therefore  higher.  In 
order  to  determine  the  impact  of  the  new  source  description 
on  the  present  sensitivity  calculation,  baseline  calculation 
No.  5 in  Table  3.2  was  repeated.  The  new  source  description 
gave  an  RDP,  \p  (°°)  , of  1.92  cm3  which  is  42  percent  higher  than 
the  old  value  of  1.35  m3. 

Results  of  the  new  sensitivity  calculations  are  pre- 
sented in  Table  3.3.  The  parameters  <p  , P , P , P , P and 

0 c ill  0 

f (or  p ) were  held  constant  at  the  values  presented  for  cal- 
w o 

culation  No.  5 in  Table  3.2.  The  near  source  seismic  velocity, 

c,  the  shear  modulus,  y,  and  the  unconfined  yield  strength,  Y , 

o 

were  varied  systematically.  Calculation  No.  13  in  Table  3.3 
is  identical  to  calculation  No.  5 in  Table  3.2  except  that  the 
new  source  description  was  used  in  Table  3.3.  The  change  in 


t 
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teleseismic  magnitude,  Am,  is  given  by  Eq.  (3.4).  Note  that 
calculation  No.  17  has  been  selected  as  the  reference  in 
Table  3.3. 

Referring  to  Table  3.3  (calculations  13,  18,  and  19) 

we  see  'that  increasing  the  unconfined  yield  strength,  Y , 

o 

decreases  the  seismic  magnitude.  This  result  is  consistent 
with  earlier  calculations  which  showed  that  increasing 
and  decreasing  Pm  both  cause  the  seismic  magnitude  to  be 
lowered.  In  general,  if  any  parameter  is  varied  such  that 
the  strength  of  the  material  is  enhanced,  the  seismic  magni- 
tude will  be  reduced.  Calculations  13,  16,  and  17  indicate 
that  increasing  the  shear  modulus  reduces  ground  motion 
coupling.  Finally,  calculations  13,  14,  and  15  show  that 
increasing  the  sound  speed  also  decreases  seismic  magnitude. 

3.3  PREDICTION  OF  TELESEISMIC  GROUND  MOTION  FROM  THE  PAHUTE 

MESA  EXPLOSION,  MAST 

Reported  here  is  a calculation,  using  the  SKIPPER  code, 
of  the  near  source  seismic  coupling  of  Mast,  a nuclear  event 
detonated  recently  in  the  rhyolite  lava  of  Pahute  Mesa.  We 
were  provided  with  a CEP  document  that  summarized  data  on  the 
material  properties  for  Mast. 

Additional  material  properties  data  were  deemed  neces- 
sary, however,  in  order  to  define  an  equation  of  state.  In 
particular,  no  information  was  available  on  shear  modulus,  bulk 
modulus,  or  failure  surface.  Since  hydrostatic  load-unload 
data  were  also  not  available,  we  chose  to  use  the  TAMEOS  fast 
table  lookup  mixture  equation  of  state  scheme.  Recognizing 
that  the  given  sound  speed  was  quite  high,  and  therefore  both 
the  bulk  modulus  and  shear  modulus  must  be  large,  a high  bulk 
modulus  was  chosen  for  the  rock  component.  The  TAMEOS  scheme 
assumes  pressure  equilibrium  between  the  rock  and  water  com- 
ponents and  the  mixture.  The  bulk  modulus  of  the  mixture,  K(^, 
is  related  to  the  bulk  modulii  of  water  and  rock  (Kw  and  KR 
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respectively)  as  follows: 


K 


w 


w 


(3.5) 


where  X and  Xn  are  the  volume  fractions  of  water  and  rock 
respectively.  For  the  given  amount  of  water  for  Mast,  at 
zero  pressure,  the  bulk  modulus  of  the  mixture  can  be  at  most 
204  kbars,  even  for  very  large  Kn,  according  to  the  mixture 
assumptions.  At  the  overburden  pressure  of  200  bars  computed 
from  DOB  and  effective  overburden  density,  the  TAMEOS  scheme 
gave  207.2  kbars  for  . 

Data  were  available  for  Halfbeak,  a shot  detonated  pre- 
viously in  a similar  rhyolite  medium,  which  showed  a Poisson's 
ratio  of  approximately  0.18.  Since  the  Half beak  and  Mast 
working  points  are  very  similar,  this  Poisson's  ratio  was 
chosen  for  Mast.  The  resulting  shear  modulus  of  168.8  kbars 
was  consistent  with  both  and  the  given  sound  speed. 

In  the  absence  of  any  data,  and  since  rhyolite  is  a 
relatively  hard,  competent  rock,  a failure  surface  identical 
to  one  used  successfully  for  Piledriver,  a shot  in  granite, 
was  chosen  for  Mast.  This  surface  was  parabolic  in  shape  and 
is  given  by 


p 

to 

1 

(i  - e\ 

p 

e / 

m 

' m r 

V m ' 

Y = Y ( 1 - §-) 
m \ em  ' 


(3.6) 


where  e^,  the  melt  energy,  is  2 x 101 0 ergs/gm,  P^  = 9 kbar,  and 
Yf  = 8 kbars.  This  surface  corresponded  to  laboratory  data 
for  granite.  In  order  to  include  macroscopic  effects  of  in 
situ  rocks,  (fracturing,  pore  pressure,  etc.),  whenever  a 
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rock  clement  exceeded  a pressure  PT  of  0.5  kbars,  the  yield 

Li 

strenqth  Y was  based  on  PT  rather  than  P.  In  spherical 

Li 

geometry,  the  von  Mises  yield  surface  degenerates  into  a 

2 

single  expression  Is  I < -r-  Y,  where  S is  the  deviatoric 

io—3  10 

stress  in  the  radial  direction.  For  any  rock  element  in 


which  pressure  had  exceeded  P , the  deviatoric  stress  was 

Li 


allowed  to  relax  toward  a lower  yield  strength  Y of  0.3 

K 


kbar  according  to  the  prescription 


K.  -!yh) 


(3.7) 


where  At  is  the  time  step  of  the  calculation  and  6 is  a 

_ _ . .... 


constant,  chosen  to  be  20  * 10  sec.  This  relaxation  is  con- 
ceptually similar  to  a Maxwell  solid  and  introduces  visco- 
elasticity into  the  material  model  for  high  stresses. 


Using  the  above  material  modeling,  a SKIPPER  calcula- 
tion was  made  for  a nominal  device  yield  of  20  tons.  Figure 
3.3  shows  peak  stress  versus  range  and  Fig.  3.4  peak  velocity 
versus  range  for  this  calculation.  In  Fig.  3.4,  a discon- 
tinuity in  peak  velocity  is  noted  at  a range  of  approximately 
18  meters,  corresponding  to  the  boundary  between  a high  stress 
viscoelastic  region,  where  relaxation  of  the  deviatoric 
stresses  is  allowed,  and  a region  where  relaxation  is  not  per- 


mitted. 


In  order  to  determine  the  seismic  coupling,  the  calcula- 
tion was  continued  until  the  response  of  the  rock-water  mixture 
became  elastic.  Displacement,  velocity,  radial  stress,  and  re- 
duced displacement  potential  (RDP)  were  monitored  at  a station 
in  the  rhyolite  located  132  meters  from  the  working  point. 
Figures  3.5  - 3.8  show  plots  of  these  quantities  versus  time. 
The  steady  state  value  of  RDP,  defined  as  nx2  where  n is 


the  static  displacement  at  position  X , was  computed  at  many 

s , 0 

stations,  all  giving  6.9  * 10  cm  . Thus,  the  material  response 


is  in  fact  elastic.  The  RDP  ip(°°),  is  a measure  of  seismic 
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Figure  3.6.  Tims  history  of  velocity  monitored  at  same 
station  as  in  Figure  3.5. 
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Figure  3. 
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. Time  history  of  stress  monitored  at  same  station 
as  in  Figure  3.5. 
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Figure  3.8.  The  reduced  displacement  potential  monitored  at 
same  station  as  in  Figure  3.5. 


coupling  and  must  be  scaled  to  the  true  device  yield  by 
multiplying  by  the  ratio  of  device  yield  to  20  tons.  The 
cavity  radius  for  a yield  of  20  tons  was  computed  to  be 
3.42  meters. 

With  the  equivalent  elastic  source  calculated  for 
Mast  we  are  now  in  a position  to  begin  computation  of  syn- 
thetic seismograms  at  specified  receiver  locations,  which 
were  recently  furnished  by  the  Project  Officer. 


3.4  CODE  DEVELOPMENT 

A new  calculational  tool  has  been  developed  to  facili- 
tate those  areas  of  the  parameter  study  which  involve  phenom- 
ena not  amenable  to  a spherically  symmetric  description.  The 
inclusion  of  free  surface  generated  interference  effects  and 
water  table  effects  both  require  a two-dimensional  description 
at  late  times.  A LINK  subroutine  has  been  developed  which 
permits  the  early  time  part  of  a given  calculation  to  be  per- 
formed with  SKIPPER,  a spherically  symmetric  stress  wave  code. 
The  SKIPPER  variables  are  written  on  tape  at  a time  before 
the  spherical  shock  reaches  a free  surface  or  submerged  layer. 
Then,  the  LINK  subroutine  reads  the  tape  and  transforms  the 
variables  to  a two-dimensional  cylindrically  symmetric  mesh 
in  the  CRAM  code.  The  remainder  of  the  calculation  is  per- 
formed with  CRAM. 

A new  crack  model  was  developed  and  checked  out  in 
the  SKIPPER  code.  The  new  model  includes  a variable  yield 
surface  above  which  all  tensions  are  forbidden  by  allowing 
cracks  to  open  in  the  appropriate  direction. 

Finally,  considerable  improvement  was  made  to  our  capa- 
bility for  generating  synthetic  seismograms.  Most  important, 
the  more  accurate  generalized  ray  theory  approach  to  computing 
the  amplitudes  of  seismic  waves  traveling  through  the  earth's 
upper  mantle  is  being  tested  and  integrated  into  our  library 
of  computer  codes. 
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IV.  CONCLUSIONS 


Considerable  progress  was  made  during  the  first  three 
months  of  this  project  on  several  of  the  tasks  specified  in 
the  work  statement  for  this  contract. 

A one-dimensional  parameter  study  which  identifies  the 
dependence  of  teleseismic  magnitude  on  near  source  material 
properties  was  carried  out.  The  major  results  of  the  material 
parameter  sensitivity  study  may  be  summarized  as  follows: 

1.  Increasing  the  air-filled  porosity  greatly 
reduces  seismic  coupling. 

2.  If  any  parameter  describing  the  yield  sur- 
face is  varied  such  that  the  material  strength 
is  reduced,  seismic  coupling  may  be  substan- 
tially enhanced. 

3.  Seismic  coupling  is  relatively  insensitive 
to  water  content;  a slight  decoupling  is 
observed  with  increasing  water  content. 

4.  Increasing  the  overburden  pressure  sub- 
stantially reduces  seismic  coupling. 

The  near  field  coupling  and  the  equivalent  elastic 
source  were  computed  for  the  recent  underground  explosion, 
Mast,  at  Pahute  Mesa.  The  next  step  is  to  generate  synthetic 
seismograms  for  this  event  at  recently  specified  receiver 
locations.  The  enhanced  computational  capabilities  for 
treating  realistic  earth  structures  will  be  exercised  in 
this  experiment. 
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